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ABSTRACT 

We present new results on the evolution of the mass function of the globular cluster 
system of the Milky Way, taking the effect of residual gas expulsion into account. 
We assume that gas embedded star clusters start with a power-law mass function 
with slope f3 — 2, similar to what is observed for the Galactic open clusters and 
young, massive star clusters in interacting galaxies. The dissolution of the clusters is 
then studied under the combined influence of residual gas expulsion driven by energy 
feedback from massive stars, stellar mass-loss, two-body relaxation and an external 
tidal field. The influence of residual gas expulsion is studied by applying results from 
a large grid of TV-body simulations computed by Baumgardt & Kroupa (2007) . 

In our model, star clusters with masses less than IO^A/q lose their residual gas 
on timescales much shorter than their crossing time and residual gas expulsion is the 
main dissolution mechanism for star clusters, destroying about 95% of all clusters 
within a few 10s of Myr. We find that in this case the final mass function of globular 
clusters is established mainly by the gas expulsion and therefore nearly independent 
of the strength of the external tidal field, and that a power-law mass function for the 
gas embedded star clusters is turned into a present-day log-normal one, verifying the 
theory proposed bv lKroupa fc Boilvl (|2002l ). Our model provides a natural explanation 
for the observed (near-)universality of the peak of the globular cluster mass function 
within a galaxy and among different galaxies. Our simulations also show that globular 
clusters must have started a factor of a few more concentrated than as we see them 
today. 

Another consequence of residual gas expulsion and the associated strong infant 
mortality of star clusters is that the Galactic halo stars come from dissolved star 
clusters. Since field halo stars would come mainly from low-mass, short-lived clusters, 
our model would provide an explanation for the observed abundance variations of light 
elements among globular cluster stars and the absence of such variations among the 
halo field stars. Furthermore, our modelling suggests a natural tendency of > IG^Mq 
gas clouds to retain their residual gas despite multiple supernova events, possibly 
explaining the complex stellar populations observed in the most massive globular 
clusters. 

Key viTords: globular clusters: general - Galaxy: formation - Galaxy: halo 



1 INTRODUCTION 

Globular clusters are among the oldest components of galax- 
ies, having formed within a few hundred Myr after the 
Big Bang (|Chabover et all 1 19981 : IVandenBerg et all I2OO2I: 
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iGratton et all 120031 '). Observations of globular cluster sys- 
tems in galaxies can therefore bring important insights 
about star formation in the early universe and the formation 
and early evolution of galaxies. One of the most remarkable 
properties of globular cluster systems is their mass distri- 
bution. Observations of the Milky Way and other nearby 
galaxies show that globular clusters follow a bell shaped 
distribution in luminosity with an average magnitude of 
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-7.3 and dispersion ay = 1.2 (|Harri j [T991I ') . For 
a mass-to-light ratio of M/Lv ~ 1.5, this corresponds to a 
characteristic mass of 1.1 ■ 10^ Mq. The peak of the glob- 
ular cluster luminosity function (GCLF) app ears to be re- 
markably similar between differen t galax ies llHarrid 199ll: 
Seeker! 1 19921 : iKundu fc Whitmord l2001al lbl: iNantais et alj 



200d}_ and also at different r adii within individual galax- 
ies j Kavelaars fc B[ancs"l997|; Iff arris. Harris fc McLaughlinI 
1 19981 : iTamura et al. 2006. ) , although evidence for a fainter 
peak of t he GCLF in dwarf galaxies has recent ly been 
reported (| Jordan et al.ll2006l : Iron den Bergiil I2OO6I ). These 
findings pose strong constraints on any theory of globular 
cluster formation and their later dynamical evolution. 

Their bell-shaped luminosity function sets globular 
clusters apart from young, massive star clusters in star- 
burst and interacting galaxies ('Whitmore & Schweizer''l99d|; 
Whitmo rc et al .1999; Fall fc Zhang 1999 ; Larsen 20oi; 
Whitmore et al.1 l2002i : Ide Griis et al] l2003h and the open 



clusters of the Milky Way and other nearby spiral galaxie s 
(|Fuente Marcos fc Fuente Marcos! [20041 : iGieles et aUl2006! ). 
which generally follow a power-law distribution over lumi- 
nosities down to the faintest observable clusters. The ques- 
tion arises whether the luminosity function of globular clus- 
ters is of primordial origin, such that the turnover mas s 
is related to a primordial Jeans mass l|Fall fc Reeslll985! ). 
or whether globular clusters also started with a power-law 
mass function and the present-day peak is due to the quicker 
dynamical evolution and preferential destruction of the low- 
mas s clusters. 

lOkazaki fc Tosa! l|l995! ') and iBaumgardtl (|l998l 'l studied 
the evolution of the Galactic globular cluster system under 
the combined influence of two-body relaxation, mass seg- 
regation and an external tidal field and concluded that an 
initial power-law mass function dN/dMc ~ with slope 

(5 ^ 2.0 can, over the course of a Hubble time, be turned into 
a bell-shaped mass function with parameters similar to what 
is observed for globul ar clusters. Using more realistic mod- 
elling of the lifetimes, IVesperinil l|l99a |2000! ). on the other 
hand, pointed out that an initial power-law mass-function 
cannot evolve to the observed bell-shaped form, but that an 
initial log-normal shape of the mass function similar to the 
observed one is preserved by the dynamical evolution, offer- 
ing an alternative way t o exp lain the near-uniformity of the 
GCLF. ' Vesperini! (|200l| ) and IParmentier fc Gilmorel (|2005| ) 
also showed that power-law initial mass functions generally 
lead to radial gradients in the peak of the mass function. 

iFall fc Zhand (|200l h found that an initial velocity dis- 
tribution which is radially anisotropic and where the amount 
of radial anisotropy increases with galactocentric distance 
is able to destroy a radi al gradient in the sha pe of the 
mass function. However, IVesperini et al.! (|2003l ) analysed 
the M87 globular cluster system and, while they found 
that a strong radial anisotropy increasing with galactocen- 
tric distance would be able to reproduce the observed con- 
stancy of the peak of the mass function at all radii, they 
also showed that such strongly radially anisotropic distribu- 
tions are inconsistent with the observed kinematics of the 
M87 globular clusters, calling into question whether glob- 
ular clusters really formed w ith power-law mass functions. 
IParmentier fc Gilmorel l|2007! ) finally studied the infiuence 
of residual gas expulsion and found that a system of pro- 
toglobular clouds with a mean mass of Mci = IO'^Mq results 



in a Gaussian cluster mass function with the appropriate 
turnover at the end of the gas expulsion phase. That shape 
is then preserved by the subsequent dynamical evolution, 
thus satisfying observational constraints. 

While the above mentioned papers considered the 
infiuence of two-body relaxation, external tidal fields 
and dynamical friction on the dissolution o f star clus- 
ters, with the except i on of iKroupa fc Boilvl (|2002l ) and 
IParmentier fc Gilmorel (|2007l ). all papers neglected the in- 
fiuence of residual gas expulsion. However, there is ample 
observational evidence that residual gas expulsion is an im- 
portant mech a nism in the early evolution of star clusters: 
iLada fc Ladal l|2003l ) for example studied embedded star 
clusters in the solar neighbourhood and found that less than 
4% - 7% of the embedded star clusters survive the initial gas 
removal to become bound clu sters o f Plei ades age. Even for 
the surviving ones, ,Weidner et al] l|2007h found that they 
lose of order 50% of their stars due to gas expulsion. A simi- 
lar steep decline of the number of c lusters with age was also 
found bv lChandar. Fall fc Whitmor c (2006 ) for star clusters 
in the SMC (see however Gieles. Lamers fc Portegies ZwartI 
(2006) who attribu te this decline to de t ection incomplete- 
ness). In addition, iBastian fc Goodwin! (|2006l ) found that 
the luminosity profiles of young, massive star clusters in 
several dwarf galaxies can best be understood by clusters 
which have undergone a rapid removal of a significant frac- 
tion of their mass as a result of gas expulsion. Finally, the 
majority of OB stars in the Milky Way and other nearby 
galax ies are found in unbound associations (|Pellerin et al.l 
[2OO6I), which, together with the assumption that stars form 
in clusters, again points to the rapid dissolution of most 
clusters within 10 Myr. Since two-body relaxation or tidal 
shocks dissolve clusters only relatively slowly, such observa- 
tions point to residual gas expulsion as the main process for 
th e dissolu t ion of star clus ters, as summarised and stre ssed 
bv lKroupal (|2005l ') (see also lBaumgardt fc Kroupal (120071 ') for 
a summary of analytic and numerical studies on residual gas 
expulsion) . 

The present paper is the eighth paper in a se- 
ries which explores the infiuence of residual gas expul- 
sion on the dynamical evolution of star clusters and 



star cluster svstems (iKroupa. Petr fc McCaughrean 19991: 
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of A'^-body simulations studying residual gas expulsion and 
the subsequent reaction of the star cluster to the drop in 
cluster potential. We varied the star formation efficiency, 
gas expulsion timescale and strength of the external tidal 
field and obtained a three-dimensional grid of models which 
can be used to predict the evolution of individual star clus- 
ters or whole star cluster systems by interpolating between 
our runs. Here, we apply these results to the globular cluster 
system of the Milky Way in order to derive constraints on 
the initial properties o f the Galactic globular c luster system. 
In a companion paper (|Parmentier et al1l2008l ) , we will focus 
on the shape of the mass function of young clusters which 
have survived gas expulsion, and study how this shape re- 
sponds to variations in the local SFE. 

In t his contribution, we ad dress two key findings de- 
duced bv lKroupa fc Boilv! (|2002! ). namely that residual gas 
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expulsion (i) re-shapes the cluster MF within ^ 100 Myr and 
(ii) that the Galactic population II stellar halo naturally re- 
sults from infant cluster mortality, thereby being physically 
related to the globular cluster population in terms of stem- 
ming from th e low-mass clusters formed at the same time 
(iLarsenlbOol ). These results were arrived at using analyti- 
cal work. Here we return to these with the help of the cur- 
rently largest existing A^-bod y library of clust er evolution 
under residual gas expulsion IjBaumgardt fc Krou pa 200X)- 
The paper is organised as follows: In Sec. [2] we describe our 
model for the initial cluster distribution and our assumption 
for the various destruction mechanisms. Sec. [3] describes the 
results and in Sec.|4]we draw our conclusions. 



a chosen global anisotropy /3„ of the cluster system, defined 
by 



/3. 



(2) 



where ve, and Vr are the two tangential velocities and 
the radial velocity of each cluster and the sums run over all 
clusters. Most simulations were done with isotropic — 0) 
or mildly radially anisotropic = 0.5) velocity dispersions. 
After setting up the clusters, we calculated the peri- and 
apocenter distances and the orbital period for each cluster 
in order to be able to estimate the influence of the various 
destruction mechanisms. 



2 THE MODELS 



2.2 Destruction mechanisms 



2.1 Initial cluster properties 

We assume that pre-cluster molecular cloud cores are dis- 
tributed with a power-law mass function dN/dMci ~ 
M^f'^' between lower and upper mass limits of Mlow = 
10^ M© and Mup = 10* Mq. Including cores with masses be- 
low IQPMq would not change our results since the clusters 
formed out of such low-mass cores are destroyed by dynami- 
cal evolution. Most simulations were made with a power-law 
index j3ci = 2.0, similar to the observed slopes for young star 
clusters in nearby and starburst galaxies. 

The star formation efficiencies e are assumed to follow a 
Gaussian distribution with a mean of 25% and a dispersion 
of 5%. Such a distribution is in agreement with observed 
star formation effic iencies which generally have e 40% 
iLada fc LadalbOOd ). No correlation of the SFE with cloud 
core mass is assumed and the gas embedded cluster masses 
are calculated according to Med ~ ^Mqi such that /3 = f3ci ■ 
We assume that star clusters are in virial equilibrium prior 
to gas expulsion. Assuming an initially cold or hot veloc- 
ity distribution for the stars would change t he impact of 
;as expulsion and the survi v al lim its (see e.g. IVerschuerenI 
1990l ): [Goodwin fc Bastiaiil (|2006l )'l. However, dynamically 
cold systems can probably be excluded for the majority of 
clusters since the strong impact that gas expulsion has on 
cluster systems would be hard to explain with them, at least 
as long as typical SFEs are of order 30%. The pre-gas ex- 
pulsion cluster radii are assumed to follow a Gaussian dis- 
tribution with width log(cri{/pc) = 0.2 and various means 
given by log(r/,/pc) = log(rh„/pc) -I- fc,. log(7?G/kpc), i.e. 
our distributions are allowed to change with Galactocentric 
distance (note that log = logj^y). 

The number density distribution of clusters in the Milky 
Way before gas expulsion is given by: 
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where Rg is the Galactocentric distance and Rcore the core 
radius of the cluster distribution, which was set to Rcore = 
1 kpc. Most simulations were done with either aa = 3.5 or 
aa = 4.5. 

We treat the Milky Way as a spherical system with 
constant rotation velocity Vc ~ 220 km/s and distribute 
the clusters spherically symmetric according to the chosen 
density profile. Cluster velocities were assigned according to 



2.2.1 Gas expulsion 

Gas expulsion was m odelled by interpolatin g betw een the 
grid of runs made bv iBaumgardt fc Kroupal |20o3). In or- 
der to use their simulation grid, the ratio of the gas ex- 
pulsion time scale to the cluster's crossing time, tm /tcross, 
the ratio of the cluster's half-mass radius to its tidal radius, 
Vh/rt, and the star formation efficiency, e have to specified. 
The star formation efficiency, e, and half-mass radius, r^, 
were already chosen when the cluster system was set up. 
The crossing time, tcross, is calculated from the pre-gas ex- 
pulsion cloud core mass Mci, and virial radius, rv, according 
to 



tCr 



2.82 
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(3) 



We assume r„ = l.SOvh, which is the relation for 
a Plummer model. If we neglect contributions coming 
from the ellipticity of the cluster orbit, the tidal radius 
of a cluster with mass Med, moving through an isother- 
mal Galactic potential with per icenter distance, Rp, is 
l|lnnanen. Harris fc Webbinkll 198^ 1 
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The model for the gas expulsion timescale tm will be de- 
scribed in Sec. 13.21 below. 

Once values for all three parameters are specified, the 
fate of each cluster was c alculated by linearly interpola ting 
between the grid points of lBaumgardt fc Kroupal (|2007l ). In- 
terpolation was done by using the 8 grid points surrounding 
the position of each cluster. Linear interpolation between 
these points in one coordinate, for example e, creates 4 data 
points with the same SFE as the cluster. A second linear 
interpolation between these 4 points in another coordinate 
reduces them to two data points and a final interpolation 
between these in the third coordinate gives a prediction for 
the surviving mass and final half-mass radius of the cluster. 
If 5 or more of the 8 grid points surrounding the position 
of a cluster correspond to dissolved clusters, the cluster was 
also assumed to be dissolved due to gas expulsion. The final 
cluster radius was determined similarly, except that interpo- 
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lation was done only if both points correspond to surviving 
clusters. 

For clusters for whic h the parameter v alues are outside 
the range considered bv iBaumgardt fc Kroupa (20o3), the 
above procedure is not applicable. In this case we assumed 
that clusters with e < 0.05 or r^/rt > 0.2 will not survive, 
as ind icated by the simulations of IBaumgardt fc Kroupal 
(|2007l l, while clusters with r^/rt < 0.01 or tm /tcross > 10.0 
were assumed to follow the same evolution as clusters with 
rh/n — 0.01 and tm /tcross ~ 10.0 respectively. 

2.2.2 Stellar evolution 

Stellar evolution reduces the masses of star clusters 
by about 30% over the course of a Hubble time 
jBaumgardt fc Makinoll2003l ). Most of this mass loss hap- 
pens within the first 100 Myr after cluster formation. We 
therefore applied stellar evolution mass loss after gas ex- 
pulsion and before the other mechanisms. Mass loss due to 
stellar evolution also causes the clusters to expand. We as- 
sumed that the expansion is happening adiabatically since 
the timescale for mass loss due to post gas-expulsion stellar 
evolution is much longer than the crossing time of star clus- 
ters, and increased the cluster radii by a factor Vhf/rhi ~ 
1/0.7 = 1.43 (.Hills 198Q ). 

2.2.3 Two-body relaxation and the Galactic tidal field 

The effects of two-body relaxation and a spherical exter- 
nal tidal field were modelle d according to the results of 
IBaumgardt fc Makinol (120031 ). who performed simulations 
of multi-mass clusters moving throu gh spherically symmet- 
ric, iso thermal galaxies. According to IBaumgardt fc Makinol 
(|2003ll , the time until disruption, or lifetime, of a star clus- 
ter moving through an external galaxy with circular velocity 
Vc on an orbit with pericenter distance Rp and eccentricity 
e is given by 



toisR. 



N 



Rp ( Vc 



ln{0.02N) J [kpc] \^220km/s 



(1 + e). (5) 



Here A'^ is the number of cluster stars left after gas expul- 
sion, which can be calculated from the cluster mass after 
gas expulsion and the mean mass of the c luster stars as 
N = Mo/<m>. A standard iKroupal l|200ll ) IMF between 
mass limits of 0.1 and 15 Mq has < m >= 0.547 Mq. x 
and k are constants describing the dissolution process and 
are g iven by a; = 0.75 and k = 1.91 (|Baumgardt fc Makinol 
I2OO3I ). An exponent x flatter than unity is also i ndicated by 
the observations of iBoutloukos fc LamersI (|2003h . 



2.2.4 Disc shocks 

Since the simulations of IBaumgardt fc Makinol (|2003h did 
not take into account the effects of passages through galac- 
tic discs, we have t o add these separately to our simu - 
lati ons. According to [Ostriker. Spitzer. "fc^ ChevaUeij (|l972l l 



and iBinnev fc Tremainel ( 19871 ). the time it takes for stars 
to increase their energy by an amount equal to their typical 
energy is given by 



tshock 



0-' Vl 



(6) 



Here is the azimuthal period of the orbit, taken to be 
equal to the orbital period, a is the internal velocity disper- 
sion of the cluster, V± is the velocity of the cluster relative 
to the disc at the time of the passage, z"^ is the average 
square-z expansion of the stars in the cluster and Qz is the 
vertical component of the gravitational field of the galaxy at 
the point where the cluster crosses the disc, and can be 
calculated from the clusters half-mass radius and mass ac- 
cording to = O.AGMo/rh and z^ = |r^. For exponential 
discs, flz can be approxirn ated by = 27vG'£oexp{R/ Rd) 
(|Binnev fc Tremainel [19871 ). where So = 750.0AfQ/pc is 
the surface density of the Galactic disc and Rd — 3.5 kpc 
its scale length. In order to obtain t shock, we calculated Vj_ 
and for pericenter and apocenter passages separately and 
summed up their contributions. The dissolution time for a 
cluster evolving under both relaxation and disc shocking was 
assumed to be 



1 
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We assumed that the mass is lost linearly over the life- 
time of a cluster, so the mass remaining at a time t < toiaT 
is given by 



Mc{t) = G.7 Mo {I ~ t/tD^aT) 



(8) 



where the factor 0.7 takes into account the mass lost from 
clusters due to stellar evolution and Mo is the cluster mass 
after gas expulsion induced loss of stars. 



2.2.5 Dynamical friction 

Massive clusters in the inner parts of galaxi es spi- 
ral in to the centres due to dynamical friction (|Spitzeij 



19691: iTremaine. Ostriker fc Spitzer! Il975l ). According to 
Binnev fc Tremainel (| 19871 ). the time to reach the galactic 
centre for a cluster of mass Mq and initial distance Rg is 
given by: 
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In our calculations, we used InA = 10. We computed for 
all clusters the friction time scale after applying gas expul- 
sion and destroyed those with tf^ic < Tnubbie. We used the 
semi-major axis of each cluster's orbit as the distance to 
the galactic centre. For the surviving clusters, their orbital 
parameters were reduced according to 



Rg f = Rg % 



tpric 



(10) 



which follows from eq. 7-28 of iBinnev fc Tremainel (|l987h . 
We assumed no change of orbital eccentricity due to dy- 
namical friction and neglected changes in the efficiency of 
dynamical friction due to the mass loss of the clusters. We 
also neglected changes in the shocking timescale due to the 
inspiral. 

We followed the evolution of the cluster system for a 
Hubble time, assumed to be TnubUe ~ 13 Gyr. Clusters 
were destroyed if either toiaT < THubbie or tpric < Tnubbie 
or if the clusters half-mass radius was larger than 33% of 
the perigalactic tidal radius after a Hubble time. 
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Figure 1. Initial (dashed lines) and surviving (solid lines) mass 
distribution of star clusters in case of no gas expulsion for clus- 
ters with Galactocentric distances Rq < 8 kpc and clusters with 
Rc > 8 kpc, co mpared to the observed distribution of Milky 
Way clusters from iHarrid l|l99(J > (points) . While the distribution 
of inner clusters is in rough agreement with the observations, too 
many low-mass clusters survive in the outer parts, due to the 
weak tidal field. 



3 RESULTS 

3.1 Runs without gas expulsion 

We first discuss the results of runs witiiout residual gas ex- 
pulsion. As a starting case, we assume a cluster distribu- 
tion with ac = 4.5, power-law mass function index /3 = 2.0 
and mean cluster radius log(rh/pc) = 0.3-1-0.2 log(-RG /kpc). 
This distribution will henceforth be referred to as the stan- 
dard case. Fig. [1] depicts the surviving mass distribution 
of clusters when we split up our sample into inner clusters 
which have Galactocentric distances Rq < 8 kpc and outer 
clusters with Rq > 8 kpc at the end of the simulation, and 
compares it wit h the sample of Milky Way globular clus- 
ters, taken from iHarrisI (| 19961 ). It can be seen that the final 
distribution for clusters inside 8 kpc (solid line) is in rough 
agreement with the observations (points), since both distri- 
butions have a maximum around 10^ Mq. This is due to effi- 
cient cluster destruction as a result of the strong tidal field in 
the inner galaxy, which reduces the number of clusters with 
masses Mc ~ Mq by several orders of magnitude. How- 
ever, the distribution of outer clusters is in contrast to the 
observations since, while the observations show a turnover 
near 10^ Mq, the simulated distribution is still rising towards 
the smallest studied cluster masses due to the weak tidal 
field in the outer parts of the Milky Way. Since the observed 
mass distribution for clusters with masses Mc > W^Mq IS 
likely to be complete, the mismatch cannot be due to in- 
completeness, but must be due to our assumptions. 

Varying the initial cluster distribution also does not 
help in reconciling this difference (Fig. [2]). For example, 
changing the cluster orbits from an isotropic distribution 
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Figure 2. Mean mass < logMc > of the model clusters compared 
to the observed distribution of Milky Way clusters for different 
initial conditions for the case of no residual gas expulsion (i.e. 
e = 1). The standard case (solid line) has ao = 4.5, logrhm = 0.3 
and /3 = 2.0 and clusters on isotropic orbits. All other runs vary 
the initial conditions; Clusters on radially anisotropic orbits with 
f3v = 0.5 (short dashed), a flatter mass spectrum of embedded 
clusters with (3 = 1.5 (dotted), clusters with larger half-mass radii 
log^'hm = 0-5 independent of Galactocentric distance (kr = 0, 
long dashed), and clusters with a flatter radial distribution in the 
galaxy qq = 3.5 (dashed-dotted). In all cases, a decrease of the 
mean mass with Galactocentric distance is predicted which is not 
observed. 



to a radially anisotropic distribution with /3„ — 0.5, or by 
choosing a flatter radial density distribution of the clusters 
in the galaxy has an almost negligible influence on the fi- 
nal mass distribution. The mismatch with the observations 
can be reduced if the initial power-law index /3 is decreased, 
since in this case a smaller number of low-mass clusters form 
initially. However, the number of high mass clusters with 
masses Mc > W'^Mq is significantly overpredicted in this 
case. In addition, observed power-law distributions have gen- 
erally /3 = 1.8 or larger. Changing the distribution of cluster 
radii to a distribution with radii < logVh >= 0.5, indepen- 
dent of Galactocentric distance, also has nearly no infiuence 
on the final mass distribution. 

All models discussed so far predict a strong change of 
average cluster mass with Galactocentric radius due to the 
fact that cluster dissolution through either relaxation or ex- 
ternal tidal shocks depends on the strength of the exter- 
nal tidal field. This variation is however not observed in ei- 
ther the Milky Way or ext ernal galaxies (jKavelaars fc HanesI 
119971 : [Jordan et al.l [ioosi ). at least not with the predicted 
strength. Hence, the observed peak cannot be caused by 
these processes but must be either primordial or due to a 
different dissolution process. 

These conclusions agree with IVesperinil (|200ll ) and 
iParmentier fc Gilrnorj ll2005h . but are in disagreement 
to iFafl fc Zhand l|20oil) and iMcLaughlin fc Falll l|200'i1) . 
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The reason for this discrepancy is that iFall fc Zhan3 
l|200ll ) studied clusters on highly radial o rbits which are 
ruled out obscrva tionally a, t l east f or M87 (jVesperini et alJ 
I2OO3I '). [McLaughlin fc Falll (|2007l ') on the other hand 
assumed dissolution times which are independent of 
the strength of the external tidal field, which is ruled 
out by simulations which s how that isolated clus- 
ters don't dissolve at all l|Baumgardt. Hut fc Heggid 
|2002| ). while for clusters in tidal fields the dissolu- 
tion time depends on the strength of the external 



tidal field ("Gne din fc Ostriker 
19971: 



[Ba umgar di I2001L 

Lamers. Gieles fc Portegies Zwart 20051 ). 



19961: Ves perini fc Heggij 



Baumgardt fc Makind l20oi 



As already explained in the Introduction, gas expulsion 
is an excellent candidate for a process which is (nearly) in- 
dependent of the strength of the external tidal field, and we 
will study its influence in the next section. 

3.2 A model for residual gas expulsion 

Neglecting the influence of the surrounding interstellar 
medium (ISM), gas from a cluster can only be lost if the 
total energy put into a gas cloud by OB stars exceeds the 
potential energy of the cloud. The potential energy of a gas 
cloud with radius r, mass Mci and star formation efficiency 
e is given by 



Eaas = k{l~e)- 



^ci 



(11) 



where the factor (1 — e) accounts for the fact that some 
gas was transformed into stars. For r = and a Plum- 
mer model, the dimensionless constant k is approximately 
fc ~ 0.4. The energy Eoas is the minimum energy which 
has to be injected into a gas clo ud in order to disperse 
it. Ac cording to the simulations bv lFreyer. Hensler fc Yorkd 
l|2006ll , who studied energy-deposition of stellar feedback 
into gas clouds, a star of mass m — SSAf© puts an energy of 
E = 6.7- erg/Myr into the ISM in the form of radiation 
and mechanical energy. Corresponding values for m = 60Mq 
(SSMq) stars are E = 18 ■ lO '^ " (E = 3.4 ■ 10^°) er g/Myr 



I Frever. Hensler fc Yorkd l2003l : iKroeger et al] |2007| ). The 

above values can be fitted with 

logio ^./erg/Myr = 50.0 + 1.72 • (log^, m/Mo - 1.55) .(12) 

Integratin g this over all stars in a cluster, assuming a canon- 
ical IMF ljKroupall200ll l between 0.1 and 120 Mq, gives for 
the total energy input of a cloud of mass Mci'- 



E 



eMc 



E,N(m)dm, 



2.5 • 10* e Mci MekmVsVMyr 



(13) 



If we compare this with the total potential energy of the 
cloud, we obtain as an estimate for the gas expulsion 
timescale tm (defined in eq. 3 of iBaumgardt fc Kroupal 
|2003)): 



TM 



Ecas/E = 7.1 ■ 10" 



1 - e Mci 
6 [Mo] 



\PC]) 



Myr. (14) 



In very massive systems, the above formula leads to gas 
expulsion timescales in excess of 3 Myr, so supernova explo- 
sions also become important in removing the gas. If we as- 
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Figure 3. Energy of a gas cloud (solid line) vs. cloud mass Mci in 
comparison with the energy input from OB stars (short dashed) 
and the energy input by supernova explosions (dotted and long 
dashed) for a gas cloud with a half-mass radius r^j = 0.5 pc and 
SFE of 25% (top panel) and the resulting gas expulsion time scale 
TM (bottom panel) used in this paper. OB stars eject enough 
energy to disperse gas clouds with masses up to 10® within 
3 Myr. More massive clouds are dispersed by supernova explo- 
sions. For the most massive clouds with Mci > IO'^Mq, the po- 
tential energy of the gas cloud can exceed the combined energy of 
all supernova explosions, meaning that some fraction of the super- 
nova ejecta will be retained. This raises the possibility of multiple 
star-formation events in these systems. The resulting gas expul- 
sion timescales for a star cluster with e = 0.25 and r/^ = 0.5 pc 
vary between 5 • 10* yrs and 3 Myrs. 

sume that a typical supernova explosion ejects around ITO^^ 
erg into the ISM, the total energy ejected by supernovae is 
given by: 

EsN = 5.025 • 10^ e fsN Mci Mekm^s'' . 



(15) 



Here fsN is the fraction of all stars with mass m > 8 M© 
and which undergo supernova explosions. For a Kroupa IMF 
from 0.1 to 120 Mq, fsN = 0.0061. Fig.|3]depicts the energy 
input of OB stars and supernova explosions into the ISM for 
gas embedded clusters with a half-mass radius of 0.5 pc and 
SFE e = 0.25. It can be seen that OB stars eject enough en- 
ergy to disperse gas clouds with masses up to 10^ Mq within 
3 Myr, i.e. before the first supernova explosions go off. The 
first supernova (assumed to have an energy of 5T0^^ erg) go- 
ing off in a cluster ejects an energy amount small compared 
to what OB stars already injected and is thus unlikely to 
have a strong impact on the cluster gas. However, the com- 
bined energy input from all supernovae is about one order 
of magnitude larger than what OB stars eject and is strong 
enough to disperse even clouds with masses up to lO"^ Mq. 

For the most massive clouds, the total energy injected 
into the ISM is however only slightly larger than the cloud 
binding energy. Some gas from supernova explosions might 
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Figure 4. Mass distribution of star clusters before gas expulsion 
(dashed lines, the "embedded cluster MF"), after gas expulsion 
(dotted lines, the "initial cluster MF") and after a Hubble time 
(solid lines), compared to the observed distribution of Milky Way 
clusters (points). Most low- mass clusters are destroyed by residual 
gas expulsion. The distribution of surviving clusters is now in 
good agreement with the observed one for both inner and outer 
clusters. 



therefore be retained in such clouds and form a second 
generation of stars enriched in heavy elements. Interest- 
ingly, in the mass range Mc > few ■ 10^ Mq multiple 
stellar pop ulations are observed i n Local Group clu sters 
like ui Cen jHilker fc Richtie3l2000l : IPiotto et all 120051 ) and 
Gl (|Mevlan et al.ll200ll ). It is also the mass range where 
a transition from globu l ar clusters to UCDs is ob served 
l|Evstigneeva et aL I l2007l : iDabringhausen et all [ioosl ). Ac- 
cording to Fig. 131 insufflcient gas removal and multiple star- 
formation events could be one reason for the differences of 
heavy clusters to ordinary GCs. 

Based on Fig. [3] we therefore choose tm from eq. [14] 
for the gas expulsion timescale if tm is smaller than 3 Myr. 
For larger gas clouds we assume tm = 3 Myr. For low-mass 
clouds, the finite time the gas needs to leave the cluster also 
becomes important. If we assume that the gas is leaving 
with the sound speed of the ISM (vs « 10 km/s), we obtain 
a lower limit for the gas removal timescale of 



TM = — = 0.098 -P^ Myr 

Vs [pc] 



(16) 



Fig. [3] depicts the resulting gas expulsion timescales for 
a gas cloud with e = 0.25 and Vh = 0.5 pc. 



3.3 Runs including residual gas expulsion 

We make the following assumptions for our standard case: 
SFEs follow a mass-independent Gaussian distribution with 
a mean of 25% and a dispersion of 5%. The effect that 
other SFE distributions have on the shape of the final clus- 
ter distribution will be discussed in more detail in a forth- 
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Figure 5. Mean masses < logMc > of the simulated clusters in 
runs with gas expulsion compared to the observed distribution of 
Milky Way clusters. The standard distribution (solid line) pre- 
dicts average masses which are in agreement with observations at 
all Galactocentric distances. Assuming radial anisotropic orbits 
(short dashed) makes little difference to the final mass function. 
Larger initial cluster radii decrease the average masses of clusters 
at all radii (long dashed) while steeper embedded mass-functions 
lead to too many low-mass clusters surviving at large radii (dot- 
dashed) and flatter embedded mass- functions overpredict the av- 
erage cluster mass at all radii (dotted). 



coming paper (jParmentier et al]|2008l ). The spatial distri- 
bution of the clusters in the galaxy has a power-law ex- 
ponent of ac ~ 4.5, the power-law mass function index 
is /3 = 2.0 and mean initial cluster radii follow a relation 
logio Th/pc = -0.1 -f 0.2 ■ log^p Ro/kpc. 

Fig. Udepicts the evolution of the mass function of clus- 
ters if we include gas removal into the runs. After gas expul- 
sion, the number of low-mass clusters is decreased by a factor 
of 10 to 100, but still continues to rise towards low-masses. 
Massive clusters are not strongly effected by gas expulsion 
due to the large ratio of ru/t cross, so their number is close 
to the initial one and the slope at the high-mass end fiattens 
only slightly. In total, only 3% of all clusters survive residual 
gas expulsion in this case. Due to the efficient destruction of 
low-mass clusters, the overall mass function after a Hubble 
time is in good agreement with the observed one for both 
inner and outer star clusters. 

The main remaining difference to the observations is 
that our runs overpredict the number of high-mass clusters. 
This could be removed in a number of ways, like assuming 
shorter gas removal times for high-mass clusters or that the 
mass function of high-mass clusters in the Milky Way was 
truncated around 10* M0 {Gieles et al, ,200^ ). due, for ex- 
ample, to star formation rates which were not high enough to 
allow the formation of more massive clusters (jWeidner et al.l 
Also, the m ass function of embe dded clusters could 
have been steeper (jWeidner et al.|[2004l ). 
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The good agreement is confirmed by Fig. [5] which shows 
the average cluster mass as a function of Galactocentric dis- 
tance. For our standard model, the average cluster mass 
decreases only slightly with Galactocentric distance, and 
stays for most distances within the limits observed for Milky 
Way globular clusters. The dispersion is also close to the 
observed one and exceeds it only at small Galactocentric 
radii. Changing the velocity dispersion of the globular clus- 
ter system from isotropic to radially anisotropic orbits with 
/3„ = 0.5 (dashed lines) leads to average cluster masses which 
are slightly larger than in the isotropic case. The differences 
are however small. Similarly, a distribution with initial slope 
P — 2.3 leads to only small changes in the final distribution. 
An embedded mass function with (3 = 2.3 might therefore 
also be compatible with the data and might in fact resolve 
the problem of too many clusters at the high-mass end noted 
in Fig. |4] Changing the distribution of cluster radii (long 
dashed line) leads to a slightly worse fit at large Galacto- 
centric distances but to a better fit in the inner parts. A 
distribution with P — 1.5 has too many high-mass clusters 
almost everywhere. The Milky Way globular cluster system 
should therefore have started with a power-law exponent in 
the range j3 « 1.8 — 2.3. Similarly, there is some room for 
variation in the other parameters since most of them lead to 
acceptable fits. 

Fig. |6] compares the distribution of cluster radii with 
the observations. Here, we assume that the projected ra- 
dius rhp is related to the three-dimensional radius vh ac- 
cording to Thp = 0.73 rh, close to the empirical relation for 
many King profiles. We also neglect mass segregation and 
assume that cluster mass follows cluster light. This should 
be a valid assumption since most globular clusters are not 
in core-collapse. 

In our standard model, the average cluster radii were 
initially increasing with Galactocentric distance according 
to < logr^/pc >= —0.1 -I- 0.2 • logi?G/kpc. Due to gas ex- 
pulsion, the clusters expand on average by a factor of 2 to 
3 and subsequent stellar evolution mass loss leads to an- 
other adiabatic expansion by 30%. It can be seen that the 
resulting distribution of cluster radii is in very good agree- 
ment with the observed distribution for both inner {Re < 8 
kpc) and outer clusters. Our adopted slope of kr — 0.2 
leads to a final slope which is clo s e to the value of 0.42 
found bv lMackev fc van den Berghl (|2005l ) in their analysis 
of the present-day Galactic globular cluster system. In this 
scenario, the half-mass radii of inner globular clusters in 
the gas embedded phase are close to the half-mass radii of 
young, embedded cl usters in the Galactic disc, ~ 0.5 pc 
jLada fc Ladal 120031 '). 

A distribution with constant logr^m/pc = 0.1 indepen- 
dent of Galactocentric distance also provides an acceptable 
fit, since a larger fraction of clusters with large radii sur- 
vive in the outer parts, leading to an increase of the average 
cluster radius with Galactocentric distance. The observed in- 
crease of half-mass radius with Galactocentric distance could 
therefore be either primordial or due to gas expulsion and 
later dynamical evolution. 

Fig.[7]finally depicts the mass density profile of the glob- 
ular cluster system as a function of Galactocentric distance. 
Gas-embedded density distributions starting from aa = 4.5 
power-law distributions lead to final distributions which are 
in good agreement with the observations (solid and dashed 
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Figure 6. Mean projected half-mass radii of the simulated clus- 
ters as a function Galactocentric distance compared to the ob- 
served distribution of Milky Way clusters in runs with gas ex- 
pulsion. The standard distribution (solid line) which has ini- 
tial cluster radii increasing with Galactocentric distance fits the 
observed distribution very well. A distribution with constant 
log^hm/pc = 0.1 independent of Galactocentric distance (dashed 
line) also provides a good fit since a larger fraction of clusters with 
large radii survive in the outer parts, leading to an increase of the 
average cluster radius with Galactocentric distance. The observed 
increase of half-mass radius with Galactocentric distance could 
therefore be either primordial or due to the dynamical evolution. 

lines). The reason is the efficient depletion of clusters at 
small Galactocentric distances, which flattens the overall 
profile. The final distribution is nearly independent of the 
initial amount of anisotropy. If clusters start with a density 
profile with power-law index aa = 3.5, similar to what is ob- 
served now, the final profile becomes too strongly fiattened 
and does not fit the observed profile. The Galactic globular 
cluster system should therefore have started more centrally 
concentrated than as we observe today and many clusters 
were lost from it in the inner parts of the Galaxy. 



3.4 Dependence of results on the assumed model 
for gas expulsion 

In order to test how these results depend on the assumed 
model for gas expulsion, we also tried two additional models 
for ga s expulsion. In one model, we follow Kroupa fc B oilvl 
(|2002l ) and assume that the gas leaves star clusters with the 
sound speed of the interstellar medium, Vs ~ 10 km/s, and 
set the gas expulsion timescale equal to tm ~ r^/vs for all 
clusters. Since observations show that star clusters are free 
of their primordial gas after one Myr (lKroupal l2005l). we set 
Tm = 0.5 Myrs in the third model. 

Fig. |S] depicts the resulting mass function of star clus- 
ters after a Hubble time for all three models. It can be seen 
that the distribution of star clusters inside Rg = 8kpc is 
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Figure 7. Final mass distribution as a function of Galactocen- 
tric distance of the simulated clusters compared to the observed 
distribution of all Milky Way clusters in runs with gas expul- 
sion. Distributions with an initial power-law profile with exponent 
aQ = 4.5 (solid and dashed lines) fit the observed distribution 
very well. Distributions with an initial power-law index oq = 3.5 
(long-dashed and dotted), which is similar to the present-day ob- 
served profile, lead to final profiles which are too flat. The Galactic 
globular cluster system must therefore have started with a more 
concentrated distribution then what is observed now. In addition, 
the orbital anisotropy seems to have only a small influence on the 
final density profile. 



always in good agreement to the observed distribution and 
does not depend much on the assumed model for gas ex- 
pulsion. The reason is the strong dissolution of star clusters 
which causes most low-mass clusters to have come from rel- 
atively high-mass, Mc ~ 10^ M©, progenitors. Since most 
of these survive gas expulsion, the final distribution is nearly 
independent of the assumed model for gas expulsion. 

In the outer parts, a factor of 5 to 10 too many low- 
mass clusters survive in case of gas expulsion with fixed 
gas expulsion timescale. But even in this case there is a 
strong depletion of star clusters by about a factor 100 com- 
pared to the initial number of clusters, so small changes 
in the assumed cluster dissolution model might bring the 
expected distribution into agreement with the observations. 
The model with constant outflow velocity on the other hand 
is in good agreement with the observations. A range of gas 
expulsion models might therefore be able to turn an initial 
power-law mass function for the gas embedded clusters into 
a Gaussian for the present-day clusters. 



3.5 The origin of the Galactic halo stars 

Using our best-fitting models, we now turn to the connec- 
tion between the Galactic halo stars and globular clusters. 



iBaumgardd (| 19981 ) and iKroupa fc Boilvl ||2002D suggested 
that if stars always form in clusters ranging from the low- 
est to the highest masses, then infant mortality and clus- 



Figure 8. Mass distribution of star clusters after a Hubble time 
for the standard model (solid lines), a model which assumes gas 
outflow with a velocity of Vs = 10 km/s (dashed lines), and a 
model which assumes a fixed gas expulsion timescale of tm = 
5 ■ 10^ yrs (dotted lines). In all cases, fewer low-mass clusters 
survive compared to a model without gas expulsion (see Fig. 1). 
The final distributions are in agreement with the observations for 
the standard case and models with fixed outfiow velocity. Models 
with fixed expulsion timescale predict a factor 10 too many low- 
mass clusters at large Galactocentric radii. 



ter dissolution will naturally produce field populations de- 
rived from the disrupted low to intermediate-mass clusters. 
A common origin of field halo stars and stars in globular 
clusters is also indi cated by simi larities in their heavy el- 
ement abundances l|Snedenll2005l ). In our standard model, 
only 5.2% of all stars initially in star clusters end up in 
present-day globular clusters. The majority of stars is lost 
from clusters either due to gas expulsion or the later dy- 
namical evolution and dissolution of the clusters. Our sim- 
ulations show that both processes contribute with approxi- 
mately equal stre ngth to stars lost from cluster s. 

According to I Freeman fc Bland-HawthornI (|2002l ). the 
total mass of the stellar halo is ~ 10^ Mq . The total mass of 
metal-poor globular clusters with [Fe/H]< — 0.8, which are 
believed to be connected to the stellar halo l|Zinnlll993l ). IS 
2-10^ M0. The initial mass in star clusters with mass Mc > 
10'' Mq is therefore about half the stellar halo mass. Taking 
into account clusters with masses 10 Mq < Mc < Mq 
would increase this estimate by another 50%. It hence seems 
possible that the entire stellar halo for med from di s solved 
star clusters, confirming t he notion of iBaumgardtl (Il998f) 
andlKroupa fc Boilvl l|2002l ). see also iParmentier fc Gilmorg 
||2007| V 

The density of field halo stars can be well fitted by a 
power-law wit h exponent a = —3. 55 ± 0.13 between 8 < 
Rg < 35 kpc l|Chiba fc Beersll2000l ). If the field halo stars 
come from dissolved, low-metallicity clusters, these clusters 
must have followed the same initial density profile. Fig. |9] 
compares the density profile of halo clusters with [Fe/H]< 
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Figure 9. Observed mass density distribution of halo clusters 
with [Fe/H] < —0.8 as a function of Galactocentric distance com- 
pared to the final distributions for clusters starting from different 
initial density profiles. The distribution starting with an initial 
power-law profile with exponent ac = 4.5 (solid line) for the 
gas embedded clusters fits the observed distribution very well. 
However, distributions with initial power-law indices in the range 
3.5 < 0(3 < 3.9, which are in better agreement with the observed 
distribution of halo stars, also provide acceptable fits to the ob- 
served cluster distribution at radii between 5 < Rq < 20 kpc. 
It therefore seems possible that halo stars and globular clusters 
share a common origin. 



0.8 with the predicted profiles, given diflFerent initial power- 
law distributions. It can be seen that density distributions 
starting from steep power-laws generally provide the best fit 
to the observations. However, density profiles starting with 
flatter initial distributions in the range 3.5 < ac < 3.8 
also provide acceptable fits to the density of halo globular 
clusters between 5 kpc < Rq < 20 kpc. It therefore seems 
very likely that field halo stars come from dissolved globular 
clusters. 

The field halo stars have a radially elongated velocity 
ellipsoid with co mponents (a'u,o'v,o'w) = (141 ±11, 106 ± 
9,94 ± 8) km/s ijChiba fc Beerd 120001 1 . corresponding to 
Pv ~ 0.5. Due to the efficient destruction of star clusters 
on radial orbits, the anisotropy of the globular cluster sys- 
tem is decreasing with time. If field halo stars come from 
dissolved star clusters, and if star clusters start with a ra- 
dial anisotropic velocity dispersion with Pv = 0.5, we predict 
a nearly isotropic distribution with /3„ = 0.01 for the surviv- 
ing star clusters after a Hubble time. The velocity dispersion 
of Galactic globular clusters can be tested if accurate proper 
motions become available, which should be possible with fu- 
ture astrometric space missions like e.g. GAIA. 



4 CONCLUSIONS 

We have followed the evolution of the Galactic globular clus- 
ter system under the infiuence of various dissolution mech- 
anisms, taking into account the effect of residual gas ex- 
pulsion. We assumed that gas embedded star clusters start 
with power-law mass functions, similar to what is observed 
for the Galactic open clusters and young, massive star clus- 
ters in interacting galaxies. The dissolution of the clusters 
was then studied under the combined influence of residual 
gas expulsion, stellar mass-loss, two-body relaxation and an 
external tidal field. The influence of gas expulsion was mod- 
eled by using a large set of A'^-body models computed by 
iBaumgardt fc Kroupal (120071 ). 

We flnd that residual gas expulsion is the main disso- 
lution mechanism for star clusters, destroying about 95% of 
them within a few 10s of Myr. It is possible to turn an ini- 
tial power-law mass function for gas embedded clusters into 
a present-day log-normal one, because clusters with masses 
less than 10^ lose their residual gas on a timescale shorter 
than their crossing times, as shown by our feedback analysis. 
In this case, the final mass function of globular clusters is 
established mainly by the gas expulsion and therefore nearly 
independent of the strength of the external tidal field, pro- 
viding a natural explanation for the observed universality 
of the peak of the globular cluster mass function within a 
galaxy and among different galaxies. Observ ational evidenc e 
for such rapid gas expulsion is discussed by iKroupal l|2005l ). 
In such a case, a characteristic mass- scale of « IO'^Mq, 
as was suggested by iFaU fc ReesI l|l985h , does not exist for 
globular clusters, and they instead form from a feature-less 
power- law m ass function, as present epoch clusters are ob- 
served to do l|Larsenl 12004 ). 

Our best-fitting model for the distribution of Galactic 
globular clusters has the following parameters: The slope 
of the mass function of gas embedded clusters is /3 = 
2.0, close to what is observed for young star clusters in 
starburst galaxies or the open clusters in the Milky Way 
dWhitmoro & Schweizer"l995"; 'Whitmore et al."'l999'; 'Larsenl 
2002; Fucntc Marcos & Fucntc Marcos 2004; Gieles et ai] 
20061 ). This slope is also similar to that found for the 
mass function of molecular-cloud cores (ISolomon et al ]|l987l : 
lEngargiola et al.1 l2003l : lRosolowskvll2005l ). Slightly steeper 
mass functions with /3 = 2.3 are also allowed, but fiatter 
mass distributions of gas embedded clusters with /? « 1.5 
lead to final distributions which are inconsistent with the 
observed mass distribution of Milky Way globular clusters. 

The current distribution of half-mass radii of globular 
clusters can be fitted with an initial log-normal distribution 
with mean log(rh/pc) = —0.1 + 0.2 log(7?G/kpc) and dis- 
persion aR = 0.2, or a distribution with log(r(i/pc) — 0.1 
and dispersion an = 0.2 independent of Galactocentric 
distance. Globular clusters should therefore have started 
with half-mass radii several times smaller than as we see 
them today. The observed dependen ce of mean half-mass 
radius on Galactocentric distance (Ivan den Berghl Il994l : 
iMackev fc van den Berghll2005^ is either primordial, point- 
ing to differences in the formation of the clusters, or coming 
from the gas expulsion and the higher survival probability 
of extended star clusters at larger Galactocentric distances. 
In terms of a unifying cluster formation theory, this latter 
solution appears more attractive. 
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The present-day radial density distribution of clusters 
can be well fitted by an initial power-law density distribu- 
tion with slope aa = 4.5 of the gas-embedded clusters. Low- 
metallicity globular clusters with [Fe/H]< —0.8 could also 
have started from flatter density distributions with slopes 
qg as low as aa = 3.5, whi ch is similar to the d ensity dis- 
tribution of halo field stars (jChiba fc Beerslbood '). 

ft seems possible that all halo field stars originate from 
dissolved star clusters. In this case halo field stars come 
mostly from low-mass clusters {Mc < 10'* Mq), which 
quickly dissolve due to either residual gas expulsion or dy- 
namical cluster evolution. In addition, gas lost by stellar 
winds from stars in low-mass clusters is likely to escape due 
to the low escape speeds. Chemical enrichment due to mas- 
sive stars within the same cluster is therefore unlikely to 
have taken place and, if the material out of which the clus- 
ters formed was well mixed, coul d explain why no abundanc e 
anomalies are seen in halo stars (jCratton et al.ll200d . [20o3 ). 

Stars in globular clusters on the other hand could 
have been enriched by material lost from more massive 
stars within the same cluster, especially if the winds from 
these stars are slow as has been recently suggested by 
iDecressin et al.l (|2007| ). Slow winds would also be ineffi- 
cient in dispersing gas clouds and therefore increase the 
gas expulsion time scale tm and hence the survival chances 
of star clusters. In order to explain the near homogeneity 
of heavy elements in globular clusters, ejecta from super- 
nova explosions should not be retained, hence the clusters 
would either have to be already gas free by the time the 
supernova explosions go off or the gas was driven out by 
these explosions. Hence, most of the enrichment of stars 
in globular clusters should take place within a few Myr, 
pointing to heavy main-sequence stars with masses between 
10 — 120 Mq as th e polluters. Thi s poin t has recently 
been stressed also by IDecressin et al.l (|2007| ). Heavy mam- 
sequence stars as sources for metal-enrichment also have the 
advantage that a standard stellar IMF would be sufficient 
and one does not have to invoke an unusual flat IMF as 
in case of poflution by AGB s tars (|D'Antona fc Caloill2004 
iPrantzos fc Charbonneill2006l '). 

Due to the efficient destruction of star clusters on ra- 
dial orbits, the velocity anisotropy of the globular cluster 
system is decreasing with time. If field halo stars come from 
dissolved star clusters, and if star clusters start with a radial 
anisotropic velocity disper sion with /3„ = 0.5 sim ilar to what 
is seen for field halo stars (|Chiba fc Beersll2000l ). we predict 
an isotropic distribution for the surviving star clusters after 
a Hubble time. 
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